Enhanced polarization in strained BaTiOa from first principles 
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ABSTRACT 

The structure, polarization, and zone-center phonons of bulk tetragonal BaTiOs under 
compressive epitaxial stress are calculated using density functional theory within the local 
density approximation. The polarization, computed using the Berry-phase formalism, in- 
creases with increasing tetragonality and is found to be enhanced by nearly 70% for the 
largest compressive misfit strain considered (-2.28%). The results are expected to be use- 
ful for the analysis of coherent epitaxial BaTiOs thin films and heterostructures grown on 
perovskite substrates having a smaller lattice constant, such as SrTiOs. 

INTRODUCTION 

With recent advances in oxide epitaxy, coherent growth of nearly perfect ultrathin per- 
ovskite films and heterostructures is now possible [|l], 0, |^. In these systems, individual 
constituent layers are under mixed mechanical boundary conditions: constrained to the 
substrate lattice constant in-plane, the films are free to relax in the normal direction. The 
strains achieved in this way can be quite significant, and can be controlled to tune structure, 
ferroelectric, and dielectric properties of these nanostructured materials |^. For example, 
the polarization could be increased over that of known bulk compounds ^ 0, improving 
the performance of ferroelectric field-effect transistors and nonvolatile memories. 

When the interfaces in lattice-matched heterostructures are between sufficiently similar 
materials, even the thinnest layers can be well described by considering the behavior of 
the bulk under the same boundary conditions. Thus, for example, many properties of 
BaTiOa/SrTiOs superlattices can be understood by considering the effects of large strains 
on bulk BaTiOs. This has led to great interest in information about highly strained bulk 
and film states in the perovskites, for example as presented previously in a phenomenological 
Landau- Devonshire framework [0, § 0. 

Here we perform first-principles calculations and obtain information about the structure, 
spontaneous polarization, and zone-center phonons of bulk tetragonal BaTiOs with misfit 
strains down to about -2.3%. While our results correspond to zero temperature, previous 
work [0] suggests that in this high-strain regime, the tetragonal phase is stable with increas- 
ing temperature and thus our analysis is relevant to the study of ultrathin BaTiOs films and 
BaTiOs/SrTiOs superlattices at room temperature |]10| . 



METHODOLOGY 



To predict the ground state structure of strained BaTiOs and compute the associated 
polarization, we use density functional theory (DFT) within the local density approxima- 



tion (LDA) ||T2|, as implemented with a plane- wave basis in the Vienna ab initio Simulations 
Package (VASP) [|l3l. Results are obtained using projector- augmented wave (PAW) po- 
tentials provided with VASP. The PAW potentials explicitly treat 10 valence electrons for 
Sr (45^4/55^), 12 for Ti {3s^3p^3d^As^), and 6 for oxygen (25^2/). The ions are steadily 
relaxed toward equilibrium until the Hellmann-Feynman forces are less than 10~^ eV/A. 
Brillouin zone integrations are performed with a Gaussian broadening of 0.1 eV during all 
relaxations. All calculations are performed with a 6x6x6 Monkhorst-Pack k-point mesh 
and a 50 Ry plane-wave cutoff. Once ground-state structures are obtained, polarizations are 
determined from first principles using the Berry-phase formalism [|15]. The force- const ant 
matrices are then obtained with a frozen-phonon approach. Specifically, a series of calcu- 
lations is performed in which each ion in turn is displaced along z by 0.1% of the lattice 
parameter and the Hellmann-Feynman forces on all the atoms computed; the phonon fre- 
quencies and eigenvectors are then acquired through diagonalization of the corresponding 
dynamical matrix. 



RESULTS 

We begin by calculating the equilibrium structural parameters of BaTiOs under zero 
stress. The observed phases of bulk BaTiOa are closely related to the cubic perovskite 
structure PmSm [|16|, which is the structure of the high-temperature phase (T>403 K). The 



simple cubic unit cell of this structure can be chosen so the Ti cation occupies the body- 
centered site lb (|,|,|) and is octahedrally coordinated by oxygen anions at 3c (|,0,|); the 
Ba cation is then at Wyckoff position la (0,0,0), and twelve- fold coordinated by oxygens. 

The rhombohedral P3ml phase is computed to be the ground state of BaTiOa, consistent 
with observations |T6[; however, it is just slightly lower in energy (only ~1 meV per unit cell) 



than the optimized tetragonal PAmm structure |T^, the observed room temperature phase 



T6|| . In the tetragonal phase, the simple cubic lattice distorts, resulting in two independent 
lattice parameters a (in-plane) and c (normal). The lower symmetry also permits relative 
displacement of the atomic sublattices along z, with Wyckoff positions la (0,0,0) for Ba and 
lb (|,|,|-|-ATi) for Ti, while two Wyckoff positions are occupied by oxygens, designated 
Oi, and On: lb (|,|,Aoi) and 2c (|,0,^-|-Ao„). These displacements result in a nonzero 
polarization along [001]. 

The computed values of the structural parameters are given in the first line of Table I. 
Our calculated lattice parameters a=3.945 A and c=3.988 A are slightly smaller than those 
found experimentally (a=3.994 A and c=4.036 A [|^), the size of the underestimate being 
typical of that found in LDA studies of a wide variety of insulators and semiconductors. 
The smaller volume is known, however, to suppress the polar instability |18|, and indeed 



the computed displacements (in reduced units of c) of 0.0122, -0.0192, and -0.0124 for Axi, 
Aoi, and Ao„, respectively, are smaller than the measured values [0, which are 0.015, - 
0.023, and -0.014. The polarization of the optimized tetragonal phase (Po=24.97 /xC/cm^) is 
thus also significantly less than that computed with the experimental structural parameters 
(32.2 fiC/cm^). As a consequence, in the following we define misfit strains to be relative 
to the calculated lattice constants, and define polarization enhancement relative to the 
computed zero-stress polarization Pq. 



In coherent epitaxial thin films and [001] superlattice heterostructures, the two lattice 
vectors in the plane of the interface are constrained to match the substrate while the third 
lattice vector is free to relax. Here, we consider tetragonal BaTiOa {PAmm) under mechan- 
ical boundary conditions corresponding to coherent epitaxy of oriented phase on 
simple cubic substrates, such as SrTiOs, having lattice constants slightly smaller than that 
of stress-free bulk tetragonal BaTiOs. We perform a series of constrained relaxations in 
which the in-plane lattice constant a of bulk tetragonal BaTiOa is compressed (e.g., to the 
value of a hypothetical substrate) and the length of the normal axis and atomic positions 
are optimized. We consider misfit strains ranging from -0.28% to -2.28%. The absolute sta- 
bilities of the resulting structures require independent investigation. For example, it may be 
energetically favorable for the polarization vector and c-axis to tilt toward the [110] direction 
(this has been referred to as the r phase in previous work lowering the point symmetry 
to monoclinic. We have, however, verified that for misfit strains of -0.76% and larger, low- 
ering the symmetry to monoclinic does not result in a lower energy, and we can therefore 
tentatively conclude that the PAmm structure is relevant under the boundary conditions 
considered in this work. 

The optimized structural parameters of strained c-phase BaTiOa are summarized in Table 
I. A tendency to conserve volume is manifested as an increase in c with increasing magnitude 
of the misfit strain, through this increase does not result in a full compensation of the in- 
plane compression and the volume actually decreases moderately. The atomic displacements 
relative to the centrosymmetrically-strained cubic perovskite structure exhibit a significant, 
nearly linear, increase with increasing misfit strain. Both the average and deviation (or 
splitting) of the two Ti-Oi distances (along the i-oriented Ti-0 chain) increase with misfit 
strain, from 1.994 ± 0.125 A at zero misfit to 2.066 ± 0.233 A at -2.28%. Thus, the smaller 
Ti-Oi distance actually decreases slightly, from 1.87 to 1.83 A, while the larger distance 
increases by 8% over the same range. The shortest Ba-On distance also declines modestly 
from 2.77 to 2.75 A. The shortest Ba-Oj distance, however, decreases more rapidly with 
increasing strain and, for the largest strain considered (2.28%), this separation is smaller 
than the minimum Ba-On contact. 

These structural trends are accompanied by a significant change in polarization Pg. 
Specifically, the increased atomic displacements associated with the increase of tetragonahty 
(c/a) tend to increase the polarization; for 2.28% strain, the polarization increases by nearly 
a factor of two. These results are consistent with a previous first-principles study, in which 
a strong dependence of the polarization on the c parameter was observed for BaTiOs [0. 
We also report the increase in energy per unit cell, AE, which reflects the energy cost of 
nonzero misfit strain. It varies quadratically with misfit strain over the range shown, with 
an estimated effective elastic constant of 135 eV/unit cell. 

The computed frequencies for the z polarized zone-center phonons for each optimized 
structure are given in Table II. As there are five atoms in each primitive cell, 12 zone-center 
optical phonons are expected (plus three acoustic modes). At the zone center these are 
labeled by the irreducible representations of the 04^, point group. There are 4 Ai modes 
(Raman and IR-active), 1 Bi mode (Raman active), and 5 doubly-degenerate E modes 
(Raman and IR-active), of which the Ai and Bi modes have atomic displacements and 
polarization (if nonzero) only along z. The frequencies of the three Ai modes are 503, 252, 
and 178 cm~^, respectively. These agree quite well with a room-temperature experiment. 



Table I: Structural parameters of BaTiOs under in-plane compression, computed from first prin- 
ciples. The misfit strain (%) is measured relative to the calculated lattice parameter of the com- 
puted tetragonal phase, reported in the table as 0.0% strain. Displacements are given relative to 
the centrosymmetrically-strained cubic-perovskite structure with the Ba position fixed at (0,0,0) 
in reduced coordinates (units of c). is the Berry-phase polarization in ^uC/cm^, with Pq being 
the polarization of the equilibrium tetragonal structure. Ai? is the increase in energy in meV per 
unit cell relative to the equilibrium tetragonal structure. 
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Table II: Computed zone-center phonon frequencies of BaTiOs as a function of misfit strain (%) 
for z-polarized phonons labelled by irreducible representations Ai and Bi. 
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which reports 512, 276, and 178 cm~ for the same TO modes ||2^. An analysis of the 
corresponding eigenvectors indicates that the two highest-frequency modes involves primarily 
Ti and O motion; in the lowest the cations move opposite to the anions, with each ion in the 
primitive cell participating almost equally. The Bi mode, for which the computed frequency 
is 294 cm~^, is Raman active and involves alternating motion of the two On ions. As for the 
Ai modes, our computed frequency agrees well with the value of 308 cm~^ measured at room 
temperature 0]. We are presently unaware of previous first-principles calculations of zone- 
center phonon frequencies in tetragonal BaTiOs, although we note previous calculations 
have been performed using an empirical rigid-ion model and shell-model potentials ||22[ 



In the remainder of Table II, we report the effect of in-plane compression on the frequen- 
cies of these four modes. We find that the two highest Ai modes increase with increasing 
misfit strain, while the frequency of the lowest-lying Ai mode is nearly independent of misfit 
strain down to about -2.3%. In contrast, the computed Raman-active Bi mode frequency is 
found to decrease with increasing misfit strain. 



DISCUSSION 



The trends witnessed here in the computed properties of strained BaTiOa are quite 
general and reflect the underlying physics of ferroelectric perovskites under the mixed- 
mechanical boundary conditions associated with epitaxy. In-plane compression resulting 
from epitaxial stress favors expansion of c normal to the substrate and, by a mechanism 
analogous to that responsible for the enhancement of the ferroelectric instability with in- 
creasing volume, the polar distortion and associated polarization along z is expected to 
increase, as observed. 

While larger polarization enhancements may in principle be achieved through larger mis- 
fit strains, the role of misfit strain in materials design is limited by a rapidly increasing strain 
energy IS.E. Above a critical thickness, misfit dislocations will become energetically favor- 
able, and the film or heterostructure will relax to its bulk lattice constant, relieving stress 
and decreasing the polarization. However, coherent epitaxy of BaTiOs-based heterostruc- 
tures of useful thicknesses has already been demonstrated 0, |^ , and the present results show 
that the misfit strain of BaTiOa on a SrTiOa substrate (~2%) is large enough to achieve 
significant enhancement of the polarization. 

Our computed phonon frequencies are of interest chiefly for two reasons. First, combined 
with the computed Born effective charges and electronic dielectric constant (to be reported 
in a separate publication), they determine the phonon contribution to the static dielectric 
response eo. While clearly an important property in itself, eo should also be relevant for 
understanding the effects of changing electrical boundary conditions on films and superlat- 
tices. Second, the sensitivity to strain of the mode frequencies, particularly the intermediate 
frequency Ai mode, can provide an experimental probe of the strain state of a BaTiOs layer. 
Indeed, for the closely related material PbTiOa, mode hardening has been observed in films 
under compressive epitaxial stress [03 . 



Although our focus here centers on the effects of strain, macroscopic electric fields, in- 
terfacial effects, and finite size effects must also play a fundamental role in determining the 
properties of coherent epitaxial thin films and superlattices. However, knowledge of strain 
effects has already proved to be invaluable for analyzing experimental and first-principles 
results |Ty], and for separating out the contributions of other intrinsic as well as extrinsic 
effects. Therefore, a broader study of BaTiOa, which includes calculation of the dielectric 
tensor and related quantities, is currently in progress. A similar study of other perovskites, 
such as PbTiOs, PbZrOa, SrTiOa, KNbOs and KTaOs, would also be pertinent. 
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